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Spatially homogeneous but possibly anisotropic cosmologies have two main types of singularities: 
(1) asymptotically velocity term dominated (AVTD) — (reversing the time direction) the universe 
evolves to the singularity with fixed anisotropic collapse rates ; (2) Mixmaster — the anisotropic 
collapse rates change in a deterministicaly chaotic way. Much less is known about spatially inhomo- 
geneous universes. It has been claimed that a generic universe would evolve toward the singularity 
as a different Mixmaster universe at each spatial point. I shall discuss how to predict whether a 
cosmology has an AVTD or Mixmaster singularity and whether or not our numerical simulations 
agree with these predictions. 
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. I. INTRODUCTION 

Einsteins theory of general relativity has passed all experimental tests to date and is thus, almost certainly, the 
classical theory of the gravitational field. Yet the universality and nonlinearity of the theory have been proven (Hawk- 
ing(1970)) to imply that regular but generic initial data will (for reasonable classical matter) develop a singularity 
either to the future or the past. In the simplest examples, this singularity is hidden from us. The Friedmann- 
^sO [ Robertson- Walker homogeneous, isotropic cosmology has a spacelike big bang hidden from our view by the cosmic 
microwave background surface of last scattering. Spherically symmetric black holes have a spacelike internal singularity 
hidden within an event horizon. Spatially homogeneous but anisotropic (different collapse rates in different directions) 
cosmologies have two types of singularities. An asymptotically velocity term dominated (AVTD) singularity means 
' that the cosmology eventually behaves as a Kasner spacetime (fixed anisotropic expansion rates characterized by a 
single parameter u) (Isenberg and Moncrief (1990)). The other type of singularity is called Mixmaster (Belinskii et 
t-H ■ al (1971), Misner(1969)). Here the spatial scalar curvature creates a potential such that a bounce off the potential 
corresponds to a change of Kasner solution. The parameter u changes in a deterministic way which is generally 
believed to be chaotic. The nature of the singularities in generic spacetimes is largely unknown. Recent analytic and 
numerical results demonstrate that some spatially inhomogeneous cosmologies approach a different Kasner singularity 
at each spatial point and thus have AVTD singularities. Belinskii, Khalatnikov, and Lifshitz (1971) (BKL) long 
ago conjectured that the generic singularity is locally Mixmaster-like. While this conjecture remains controversial, it 
!h ■ provides a paradigm which can be tested. 

Iii this talk, I shall discuss the progress toward understanding the nature of the generic cosmological singularity 
(Berger, Garfinkle, and Moncrief (1997)). In particular, I shall emphasize the application of a method due to Grubisic 
and Moncrief (1993) (GM) to predict whether one should expect the singularity to be AVTD or Mixmaster-like and 
then discuss numerical studies to support or contradict such predictions. This work was done in collaboration with 
David Garfinkle (Oakland University), Vincent Moncrief (Yale University), James Isenberg (University of Oregon) 
and graduate students Eugene Strasser (Oakland University), Boro Grubisic (Yale University), and Marsha Weaver 
(University of Oregon) . 
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II. MIXMASTER DYNAMICS 



To understand the BKL conjecture, let us first discuss a vacuum Bianchi IX cosmology. According to the Bianchi 
classification of 3-dimcnsional homogeneous spaces, the spatial line element has the form 

dl 2 = A^-oV (1) 

where the Aij are constants and the a l s (generalizations of dx l ) satisfy 

da 1 = Cj k cr- ? A <r k . (2) 

For (vacuum) Bianchi IX, C* fc = e* fe , the antisymmetric symbol. The spatial dependence of the a l s when written in 
a coordinate basis yield a non-vanishing spatial scalar curvature, 3 R. To obtain a cosmological spacetime, we shall 
let the become functions of time only. For convenience, we shall require Ay to be diagonal. It is then useful to 
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represent these anisotropic scale factors e a , , and e 7 in terms of the logarithmic volume SI and orthogonal anisotropic 
shears (3± where 

a = Sl- 2/3+, 

C = ft + /3+ + V3/3_, 

7 = SI + (3+ - V3/3_. (3) 
Einsteins equations can be obtained by variation of the Hamiltonian constraint 

H = -pi + p\ + p 2 _ +V(p + ,p_, SI) (4) 
where pn, P± are conjugate to SI, (3± and 

V(/3 + ,/3-,Q) = e 4Q ~ 80 + + e 4a + 4 /3++ 4 V^/3- + ^0+4/3+ -4^3/3- _ 2e 40+4/3 + _ 2e 4n-2/3 + -2v / 3/3- _ 2e 4f2-2^++2y3/3_ _ /g\ 

The singularity occurs as SI — > — oo. Note that Eq. (4) splits naturally into a kinetic piece 

ff* = -Pn+l£+j£ (6) 
and a potential piece. Einsteins equations also require 

H = 0. (7) 
We shall now explore the method of consistent potentials (MCP) as used by GM. 
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FIG. 1. A typical Mixmaster trajectory as projected into the anisotropy plane. The rescaled variables (Moser et al (1971)) 
cause the bounces to occur at fixed locations. More than 262 bounces are shown (see Berger, Garfinkle, and Strasser (1997)). 



The case V = corresponds in fact to the Kasner solution (Bianchi Type I with C l - k = 0). Thus Eq. (4) may be 
replaced by Hk whose variation yields in the limit SI — > — oo, 

h - gfHH. (8) 

However, we require from Eq. (7) (now Hk = 0) 

Pn=P 2 ++P- (9) 
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and can replace Eq. (8) with 



(3+ = cos9\Vl\ , /3_=sin0|fi|. (10) 

If the Mixmaster model were AVTD, substitution of Eq. (10) into Eq. (5) would yield only terms which decay 
exponentially as f2 — > — oo. 

Consider the first term on the right hand side (rhs) of Eq. (5). Upon substitution with Eq. (10), we find 

e 4fi-8/3+ _^ e -4|n|(i+2c OS e)_ 

If cos 9 < —5, this term will grow exponentially. A bounce off this exponential wall then occurs leading to a new 
value for cos#. Conservation of momentum yields (Ryan (1971)) 

, 4 + 5cos6> 

cos 9 = —- (12) 

5 + 4 cos 9 ' 

which causes the first term on the rhs of Eq. (5) to exponentially decay. However, depending on the sign of sin 9, 
either the second or third term on the rhs of Eq. (5) will begin to grow exponentially. This process repeats indefinitely. 
Numerical studies (Moser et al (1973), Rugh and Jones (1990), Berger, Garfinkle, and Strasser (1997)) confirm this 
picture (see Fig. 1). This means that all terms which were neglected in the MCP are in fact negligible. 

III. GOWDY UNIVERSES ON T 3 x R 

If spatial dependence in one direction (e.g. on < 9 < 2n) is allowed in the anisotropic scale factors of a Bianchi I 
homogeneous cosmology, the T 3 x R Gowdy (1971) solution is obtained (see also Berger (1974)). The metric 

ds 2 = e-^ 2 e^ 2 (-e- 2T dT 2 + d9 2 ) + e-\e p do 2 + 2e p Q dadS + (e p Q 2 + e - p )dS 2 } (13) 

describes amplitudes P(9,t) and Q(9,t) for the + and x polarizations of gravitational waves propagating in an 
inhomogeneous background spacetime written in terms of X(6, r). This model has the extremely nice property that 
the dynamical equations for the waves obtained by variation of the (non-zero) Hamiltonian (Moncrief (1981)) 



H = jd9 (14 + \e- 2P it 2 Q ) +jd9 (\e- 2T P 2 +^ p -^Q 2 ) = H K + H v 



(14) 



decouple from the constraints. Here 7rp and ttq are the momenta conjugate to P and Q respectively. The constraints 
are first order equations for A which can be solved trivially in terms of the P and Q found as solutions to the wave 
equations obtained from Eq. (14). 

To use the MCP, we obtain the AVTD solution as r — > 00 by variation of Hk from Eq. (14). This yields 

P^v(9)t , Q^Q a (9). (15) 

Consider substitution in the separate terms of Eq. (14). The potential 

in Hk becomes 



n=n 2 e- 2P (16) 



Vi - ir 2 Q e- 2 ^ (17) 
which decays exponentially only for v > 0. On the other hand, 

V,=Q 2 e 2 ( p -^ (18) 

in Hy becomes 

U 2 = Q> 2( "- 1)T (19) 

which decays exponentially if v < 1. Thus we see that < v < 1 (at each spatial point) is consistent asr^oo with 
an AVTD singularity. 
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FIG. 2. P, Vi, and Vi vs r at a fixed value of 9. Note that the slope of P, P, T , decreases after each interaction with V2 while 
P, r goes from negative to positive after each interaction with Vi. A continuation of this graph in r would show that Vi and V2 
have permanently died off and that P continues to increase with fixed positive slope P, T < 1. 



What happens if P, T < or P, T > 1? Then either V\ or V2 will be large. A bounce off each potential will drive 
P, T into the range allowed by that potential. Alternate bounces off both potentials will eventually drive P, T into the 
required range. This process is shown at a typical spatial point in Fig. 2 which is taken from a numerical simulation 
(Bcrger and Moncrief (1993), Berger and Garfinkle (1997)). 

We note that exceptional behavior can occur at isolated spatial points where a coefficient of V\ or V2 (i.e. ttq or 
Q,e) vanishes. Since near these exceptional points V\ or V2 is small, the mechanism to drive P, T into the range (0, 1) 
takes a very long time to become effective. This means that P, T > 1 or P, T < can persist leading to the development 
of spiky features in P and Q. These are discussed in detail elsewhere (Berger and Garfinkle (1997)). 

The MCP combined with numerical simulations provide strong evidence that the Gowdy singularity is AVTD 
except, perhaps, at isolated spatial points. 



IV. (7(1) SYMMETRIC COSMOLOGIES ON T 3 x R 

We have also considered a class of cosmologies with dependence on two spatial variables u, v as well as time r. The 
model may be described by five degrees of freedom {z, x, ip, u, A} with conjugate momenta {p z ,Px,P<p,Puj,Pa}- The 
metric takes the form (Moncrief (1986), Berger and Moncrief (1993)) 

ds 2 = e' 2 * {-N 2 e- 4T dr 2 + e- 2r e A e ab dx a dx b } + e 2ip (dx 3 + p a dx a ) 2 (20) 

where {x a } — {u, v}, N — e A , x 3 is the symmetry direction and 

! fe 2z + e- 2z (l+x) 2 e 2z + e- 2z (x 2 . . 

eafc-2 { e 2z +e -2z( x 2 _ 1} g 2z + & -2z (1 _ ^2 J W 

is the conformal metric in the u-v plane. The variable u> is obtained from a canonical transformation of the twists (3 a . 
Einstein's equations are obtained from the variation of 

H = jj dudv (\p 2 z + \e^p 2 x + \ V % + \e A ^ V l - \p\ + 2p A ) + e~ 2T j> j> dudv { (e V") , ab - (e V b ) , Q A, 6 

+e A [(e- 2z ) , u x, v - (e- 2z ) , v x, u ] + 2e A e ab tp, a <p, b +±e A e-^e ab iu, a u, b } =H K + H V . (22) 

Unfortunately, the constraints are nontrivial. So far, we only have a restricted solution to them initially and have 
not yet implemented a program to enforce them during numerical evolution. (Constraint preservation is guaranteed 
analytically but not in a differenced approximation to Einsteins equations.) 

Despite the complicated appearance of Eq. (22), the MCP may still be used. Variation of Hk in Eq. (22) yields as 
t — > oo the AVTD solution 
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z^-v z (u,v)t , x^x (u,v) , tp — » —v v (u,v)t , u^>Wo{u,v) , 
A — > A (u, u) + [2 — _Pa(u, u)]t , p z -» -4^(w,v) , p x ^p° x (u,v) , 

P v -» -4tv(u,«) , , pa^Pa(u,w) . (23) 

In the AVTD limit, the Hamiltonian constraint becomes 

H K = 2v\ + ie- 4 -(p°) 2 + 2vl + \e~^ u f » 2 = 0. (24) 

The exponential terms in Hk decay only for v z > and v v > as in the Gowdy spacetime and will act to drive z, T 
and <p, T into the allowed range. We require p\ > to enforce collapse. Thus we need only consider the restriction 

H K = 2vl + 2v 2 v -l(p°) 2 = (25) 

as t — > oo . Although Hy is complicated, there are only two types of exponential dependence. Schematically, for 
Hy = § § dudv V, 

V M Ce A - 2r - 2z + e A - 2T - 2z - 4v (VLu) 2 . (26) 
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FIG. 3. The logarithm of the non-schematic V (from Eq. (22)) vs r for a polarized U(l) model. The solid lines are exponential 
best fits. 



Only the term in Eq. (22) containing spatial derivatives of lu has the behavior shown in the second term on the rhs 
of Eq. (26). The others all behave as the first term. Substitution of the AVTD solution Eq. (23) and the restriction 
given by Eq. (25) yield (again schematically) 

Vw Ce~ 2 ^ r + e 2 "^(Vw) 2 . (27) 

The special case of polarized U(l) models have w = p^ = 0. This is the only degree of freedom that may be 
permanently removed (both analytically and numerically). In this case, 

V w Ce- 2v * T . (28) 

Thus we expect polarized models to have AVTD singularities (Grubisic and Moncrief (1994)). The numerical simu- 
lations verify this prediction (Berger, Garfinkle, and Moncrief (1997), Berger and Moncrief (1997)). Fig. 3 shows the 
exponential decay of V at two typical spatial points. Behaviors of A, z, x, and if are also consistent with an AVTD 
singularity. 

For generic models, Eq. (27) predicts a growing term in V. This means that if a collapsing spacetime tries to become 
AVTD, a potential term will grow exponentially. A bounce off this potential will (among other things) change the 
sign of <p, T so that an exponential term in Hk will grow. The process should repeat indefinitely since neither v v > 
nor v v < is consistent with both terms exponentially small. 
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FIG. 4. The potential V vs r The base-10 logarithm of the (non-schematic — i.e. from Eq. (22)) second term on the rhs of 
Eq. (27) is plotted (solid line, left scale). Note the exponentially growing and decaying portions of the graph. Also plotted 
(broken line) is the combination of momenta (pa + p v +p z /2) that yields the negative of the quantity that multiplies r in the 
AVTD limit. The horizontal line indicates zero for this quantity (right scale). 

Results from a simulation (Fig. 4) appear to be consistent with this picture which itself is in accord with the BKL 
conjecture. However, there is a large caveat. The exponents in V in Eq. (26) are regulated by the kinetic part of 
the Hamiltonian constraint Eq. (25). If the value of Hk is incorrect due to numerical error, the exponent may have 
the wrong sign yielding qualitatively wrong behavior. Thus these models will not provide information on the validity 
of the BKL conjecture until enforcement of the constraint is implemented. In polarized (7(1) models, the result is 
much stronger because the constraints remain small (and converge to zero with increasing spatial resolution) for all 
t (Berger and Moncrief (1997))). 

A cleaner example of very similar behavior in one spatial dimension will be discussed elsewhere (Weaver et al 



We explore the behavior of cosmological spacetimes in their collapse to a singularity. The MCP can predict whether 
the singularity will be AVTD or Mixmaster-like. In the Gowdy case, we learn how the nonlincarities in the wave 
equations force the parameters into the range for consistent AVTD behavior. At isolated spatial points, these terms 
are absent so non-generic behavior can occur. This leads to the growth of spiky features. 

In polarized U(l) cosmologies, these spiky features do not arise. Numerical evolutions support the prediction that 
these models are AVTD. In generic U(l) models, we expect a Mixmaster-like singularity. Numerical difficulties prevent 
a definitive conclusion in this case although there is some evidence to support this claim. 

We emphasize that, to the extent we believe the numerical results, the character of the singularity is completely 
local. The predictions of the MCP assume constancy in time of the spatially dependent coefficients of the exponentially 
growing or decaying terms. Simulations support such behavior sufficiently close to the singularity. 

While the singularity occurs at (say) r = oo, we have performed numerical simulations only to some tf < oo. 
Nonetheless, we believe that we can use these simulations to support conjectures about the nature of singularities in 
these models. Since the numerical simulations support predictions based on the MCP, we argue that we understand 
the effect of all the nonlinear terms on the evolution toward the singularity. It seems reasonable to suppose that for 
Tp < t < oo nothing qualitatively different will happen. This is almost certainly true for AVTD spacetimes. An 
infinite sequence of Mixmaster bounces at different spatial points could yield eventually a very complicated spatial 
structure (e.g. Kirillov and Kochnev (1987)). It is likely, nevertheless, that the exponential potentials will still 
dominate over increasing spatial derivatives as the controlling factor in the evolution. 
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V. CONCLUSIONS 
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